Clear["Global`*"];


r0odKordp[Kordp_]:=(
	LG=(((ΔX0)/(ΔY0))/2)yzr{-(resx/resz),1};(*levo zgorne ogl.*)
	r0=LG+((-2 LG[[1]])/resx){1,-1} (LTogl+Reverse@Kordp);  (*2D*)
	{r0[[1]],yzr,r0[[2]]}
);


(*
rinc: {pozicija,smer žarka} na meji materilov
n0, n: lomna količnika materialov
Δy: debelina novga materiala
*)
rinc[rinc00_,n0_,n_,Δy_]:=(
	{r00,c00}=rinc00;
	φ0=N@ArcTan@((Norm[c00[[{1,3}]]])/(c00[[2]]));(*sin*)
	c=RotationMatrix[φ0-ArcSin[((n0 Sin@φ0)/n)],{c00,{0,1,0}}].c00;
	r=r00+Δy (c/(c[[2]]));
	{r,c}
);

z[Kordp_]:=(  (*tuki shranmo celo 3D pot.*)
	rinc0=Table[r0odKordp[Kordp],2];(*pred steklom*)
	rinc1=rinc[rinc0,nzr,nst,yst];(*pred vodo*)
	rinc2=rinc[rinc1,nst,nvo,yvo];(*pred posodico*)
	rinc2[[1,{1,3}]]
);

(*map za shranjevanje rezultatov*)
mapa="c:\\Users\\gal\\Documents\\DRUGI PROGRAMI\\Wolfram Mathematica\\saška ura\\data iz posnetkov\\"; 
imeTp="TNji.h5";
imeLZp="vogaliLZ.h5";(*levi zgornji vogal piksel*)
imeKNp="KNji.mx";(*koord. najdenih*)

imeT="TNjim.h5";
imenagibi="φji.h5";


fps=240;(*slik na sekundo za kamero*)
res={resx,resz}={1080,1920};
(*ločljivost snemanja: širina, višina: POZOR: SNEMANJA, IN NE OBREZANEGA POSNETKA*)
ΔX0=1.6060500515995875;(*kako visoka bi bila slika na razdalji 1 od kamere (skopiramo iz )*)
ΔY0=1;(*pustimo pri meru*)
LTogl={10,693};
(*Top,Left prepisan iz premiere obrezovanja*)

yzr=22*8/1000;(*zraka*)(*22 lego*)
yst=.67/100;(*stekla*)
yvo=17.9/100;(*vode*)

nzr=1.0;(*približno*)
nst=1.5;(*približno*)
nvo=1.333;

podmape=(#<>"\\")&/@Select[FileNames[All,mapa],DirectoryQ];

Kordpjiji=Import[#<>imeTp, "Dataset1"]&/@podmape;
Kordjiji=Map[z,Kordpjiji,{2}];
Do[
	Export[podmape[[i]]<>imeT,Kordjiji[[i]]],
{i,Length@podmape}]

Kordpjijiji=Import[#<>imeKNp]&/@podmape;
nagibi=ArcTan@Map[
	A/.FindFit[#,A x+B ,{A,B},x]&,
Kordpjijiji[[;;,;;,;;,{2,1}]],{2}];

Do[
Export[podmape[[i]]<>imenagibi,nagibi[[i]]],
{i,Length@podmape}]

Lah bi tut za te računu lege, ampak brezveze: pokur vsaj 10GB rama, dolg traja, pa razlike so tko minimalne, ker je vse tok skupi: krepko krepko znotrej sicer napake
Zato tega raj ne poganjat.

Table[(*vsak posnetek*)
	Table[(*vsaka slika*)
		Kordpji=Kordpjijiji[[ip,is]];
		Table[LZpjiji[[ip,is]],Length@Kordpji]+Kordpji,
	{is,Length@Kordpjijiji[[ip]]}],
{ip,Length@Kordpjijiji}];